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A  new  family  of  life  distributions,  called  the  wear-out  distributions, 
is  developed  on  the  basis  of  a  failure  rate  function,  which  is  a  con¬ 
stant  up  to  the  change-point  and  strictly  increasing  afterwards.  Prop¬ 
erties  of  these  wear-out  distributions  are  derived  and  a  Bayes  adaptive 
procedures  is  developed  for  the  estimation  of  the  change  point.  Recur¬ 
sive  formulae  are  given  for  the  determination  of  the  posterior  probability 
that  the  change  has  occurred  and  of  its  Bayes  estimator.  The  results  of 
numerical  simulations  are  given  to  illustrate  the  properties  of  the  adap¬ 
tive  procedure. 
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1.  Introduction 


Consider  a  system  having  several  components,  whose  life  distribu¬ 
tions  depend  on  the  state  of  the  system.  Each  component  is  replaced 
immediately  after  failure  (corrective  maintenance),  or  at  scheduled 
replacement  epochs  (preventive  maintenance) .  The  optimal  scheduling 
of  the  replacement  epochs  depends  on  various  economic  considerations, 
and  on  the  particular  life  distribution  o£  the  component.  This  in  turn 
depends  on  the  state  of  the  whole  system.  We  recognize  two  phases  of 
the  system.  Phase  I  is  the  phase  in  which  the  components  have  a  con¬ 
stant  failure  rate,  X  .  Phase  II  is  the  phase  in  which  the  failure 
rate  is  greater  than  X  and  increasing.  Phase  I  is  a  mature  stable 
phase  of  the  system,  while  Phase  II  is  a  wear-out  phase.  It  is  impor¬ 
tant  to  detect  early  the  transition  of  the  system  from  Phase  I  to  Phase 
II.  Optimal  replacement  scheduling  and  inventory  management  during  the 
wear-out  phase  of  the  system  are  different  from  those  appropriate  to  the 
constant  failure  rate  phase.  The  present  paper  is  devoted  to  the  study 
of  the  detection  of  the  shift  from  Phase  I  to  Phase  II,  and  estimating 
the  change  point ,  t  ,  where  this  shift  occurs . 

In  Section  2  we  derive  from  a  specified  failure  rate  model  an 
appropriate  life  distribution  for  the  present  study.  This  life  distri¬ 
bution  is  a  composition  of  an  exponential  distribution  with  a  Weibull. 

It  seems  to  be  new  in  the  sense  of  not  being  dealt  with  in  the  litera¬ 
ture  on  reliability.  We  shall  call  it  the  wear-out  life  distribution. 
Several  properties  of  this  distribution  are  presented  in  Section  2.  In 
Section  3  we  develop  a  Bayes  adaptive  procedure  of  estimating  the  change- 
point  t  ,  on  the  basis  of  a  sequence  of  the  replacement  times  of  the 
component.  Section  4  is  devoted  to  the  derivation  of  some  recursive 
formulae,  which  are  required  by  the  adaptive  procedure.  In  Section  5 
we  provide  some  simulation  results,  which  provide  information  of  the 
characteristics  of  the  detection  procedure. 


2.  The  Wear-Out  Life  Distribution 


The  wear-out  life  distribution  is  based  on  the  following  non¬ 
decreasing  failure-rate  function 


:;X,a,T)  -  - 

► 


if  OStST, 


(2.1) 


X+xMt-T^)0"1 


where  X  is  the  constant  failure-rate  in  Phase  I,  0<X<®  ;  a  is  a  shape 
parameter,  lsa<®  and  t  is  the  change  point  parameter  T+*max(0,T)  . 

If  t£0  the  system  is  already  in  Phase  II.  The  present  model  represents 
a  superposition  of  a  Weibull  failure-rate  function  (see  Mann,  Schafer 
and  Singpurwalla  [  1  ;pp  J.27] )  on  a  constant  failure-rate  function  for 
t2x+  .  Notice  also  that  the  limiting  case  of  ct-1  corresponds  to  a 
shift  at  t  from  an  exponential  life  distribution  with  MTBF  y^*l/X  to 
an  exponential  life  distribution  with  an  MTBF  ,  i.e.,  for  all 

t>t  ,  p{X2t | x2t }  "  exp{-2X(t-r) }  .  From  this  failure-rate  model  we 
obtain  a  life  distribution,  called  the  wear-out  life  distribution,  with 


a  c.d.f. 


F(x;X,a,x)  -  1-exp {-X  /  (l+n(X(t-T.)  X)dt} 

0  + 


(2.2) 


-  l-exp(-Xx-(X(x-x+)+)  }  ,  xiO 

Obviously,  F(x;X,a,r)  *  0  for  xSO  . 

The  corresponding  p.d.f.  is 

f(x;X,r,a)  -  Xe“*x(l-h*(X(x-T+)+)°  1)  . 

( 

•exp{-Xa(x-T+)“>  ,  xiO  . 

It  follows  the  conditional  p.d.f.  of  X  (the  life  length),  given 
(X2t)  ,  which  is  the  p.d.f.  of  the  residual  life  distribution  in  Phase 
II,  is 

g(x-T+;X,a)  •  Xe”^^x“T+^  (l-hx(X(x-t+))°'*1)  . 

•exp(-Xa(x-t+)a}  ,  x2t  ( 


(2.3) 


(2.4) 


Notice  that  8*1/ X  is  a  scale  parameter  of  this  distribution,  and  that 
is  a  location  parameter  of  the  truncation  type  (see  Zacks  [  3  >PP>14Cl)> 
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LIFE  LENGTH 

Figure  1.  Density  functions  of  the  wear-out  distribution 
for  X«1,  ol  =  2  and  r*0,  1  and  co  (exponential). 
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In  Fig.  1  we  present  graphs  of  the  wear-out  life  density  (2.3)  for 
a*2,  1*1  and  various  values  of  t  .  We  remark  here  that,  if  t 

is  very  large  compared  to  1/A  ,  the  wear-out  distribution  is  practi¬ 
cally  an  exponential  distribution  with  mean  p*l/A  .  This  is  the  case 
when  the  system  is  far  from  its  wear-out  phase.  We  conclude  the  present 
section  with  a  presentation  of  the  formula  for  the  moments  of  the  wear- 
out  distribution  (2.3).  Since  1/A  is  a  scale  parameter,  we  will  pre¬ 
sent  the  moment  pr(a,T)  for  the  case  of  A-l  .  If  A*1  then 

EtX^A^a.-r}  *  ~  u  (a, At)  . 

Ar  r 

The  r-th  moment  in  the  standard  case  is  for  r>0  , 


where 


p  (a,t)  *  /  xre  xdx  +  /  xre  X  ^X  (l+a(x-t)a  ^)dx 

0  T 

■  r! (l-Pos(r |t)) 

+  e“T  £  (*)rr“JM  (a) 

j-0  2  2 

“A  k 

Pos(k| A)  »  e  A  Z  ~r 

j-0  2’ 


is  the  Poisson  c.d.f.  and 

co  Cl 

M. (a)  -  /  y^e~y~y  (l+cfy0  1)dy  , 
0 

Integration  by  parts  yields,  for  every  j£l  , 

M  (a)  -  j  /  yJ_1e“(y4y  3dy  , 
j  n 


j*0 


(2.5) 


(2.6) 


(2.7) 


(2.8) 


which  is  smaller  than  r(j+l)  -  j!  for  all  J2l  .  Notice  that  Mq(o)«1 
for  all  a2l  •  The  Integral  in  (2.8)  can  be  evaluated  numerically.  In 
the  case  of  a*2  we  obtain,  for  all  jSl  , 

H,  (2)  -  J«1/4  '  jS-K-W”2* 

J  n 


0 


.  J.1'4  I 

w**  v 


(2.9) 
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where,  for  all  tz2 


BUiS) 


•  t  12 
/  ze  2  z  dz 


5^_1exp{-  \  l2}  +  (l- 1)  B(£-2;C) 


(2.10) 


B(0;£)  -  /57  (1-4 (C)) 

(2.11) 

B(1;0  ”  exp{-  j  l2}  , 

where  4(z)  Is  Che  standard  normal  c.d.f. 

The  functions  B(£;  ~  )  can  be  determined  recursively,  according  to 

(2.10)  and  (2.11).  In  Table  2.1  we  provide  the  first  four  moments,  the 
standard  deviation  cr  ,  Che  measures  of  skewness  y^  and  kurtosis,  y^  , 
of  the  wear-out  distribution  for  A»1  ,  a*2  and  several  values  of  t  . 

Table  2.1  The  first  four  afements.  standard  deviation,  and  measure 
of  skewness  and  kurtoals  of  the  wear-out  distribution,  for  X«1  ,  a-2 
and  t»0(.5)2  ,  3(1)7  . 


0 

.5 

1.0 

1.5 

2.0 

0.54564 

0.72442 

0.83285 

0.89862 

0.93851 

0.45436 

0.78694 

1.09709 

1.35098 

1.54486 

0.47769 

1.03764 

1.76118 

2.53153 

3.25967 

0.59026 

1.55837 

3.18161 

5.34890 

7.86491 

0.39577 

<5731351 

0.63518 

0.73720 

0.81490 

0.94918 

0.65373 

0.68460 

0.85060 

1.04106 

3.80641 

2.95378 

2.68424 

2.87374 

3.32183 

0.97738 

0.99168 

0.99694 

0.99887 

0.99959 

1.78732 

1.90511 

1.95897 

1.98265 

1.99279 

4.42171 

5.15969 

5.57697 

5.79571 

5.90446 

12.93668 

17.10481 

20.00519 

21.81282 

22.85481 

0.91217 

0.96005 

0.98239 

0.99242 

0.99680 

1.38130 

1.63008 

1.79281 

1.89034 

1.94467 

4.55947 

5.85946 

6.96183 

7.77159 

8.30424 
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The  exponential  distribution  is  the  limiting  case  when  t-*»  .  The 
moments  of  the  exponential  distribution  with  X«1  are  y^«l, 
y,j*6  and  y^«24  .  We  see  in  Table  2.1  that  when  t«7  the  moments 
of  the  wear-out  distribution  are  already  close  to  those  of  the  limiting 
exponential  distribution. 
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3.  Bayes  Estimation  of  the  Change  Point 

In  the  present  section  we  develop  the  Bayesian  framework  for  estim¬ 
ating  the  change  point,  r  ,  and  deciding  whether  the  transition  to  Phase 
II  has  taken  place. 

Given  a  sequence  cx» C2’ * ’ ’  °f  replacement  epochs,  where 


o<t1<t2<. 

where  tQ  =  0 


,  we  consider  the  random  variables  X.  * 


Xi  represents  the  operational  time  length  of 


—  (i-1 ,  2 , . . . )  , 

i-th 


sequentially  Installed  component.  Since  the  preventive  maintenance 
replaces  the  component  after  A  units  of  time,  X^  *  min(Y^,A)  ,  where 
*1’X2’”*  are  independent  random  variables  having  the  wear-out  distri- 


t— t^^2  (i-1, 2, . . .) 


is  the 


bution,  with  parameters  X,o  and  t 
length  of  time  left  for  the  system,  after  the  i-th  installment  of  the 
component  till  the  transition  to  Phase  II.  If  r^<0  ,  the  chang  to 
Phase  II  has  taken  place  before  ^  . 

In  the  present  study  we  assume  that  A  and  a  are  known.  The 
only  unknown  parameter  is  the  change  point,  t  .  Without  loss  of  gener¬ 
ality  we  assume  that  X-l  .  If  X*1  replace  X^  in  the  formulae  below 
by  XX1(i«l,...,n)  . 

Given  n  replacement  points  t,<t„<...<t  ,  let  X 

l  &  n  ~ 


(*>, 


(X, . X  V 


where  X.-t.-t.  .  , 
i  i  i-1 


Ji" 


and  define 


,  if  X±  <  A 


,  if  Xt  -  A 


(3.1) 


The  likelihood  function  of  t 


r(n), 


given  X(n)  , 


is 


.0-1. 


L(t;Xv“')  -  KtsO}  erp{-  E  XT}  H  (l+aJ,X7  A) 

i-1  1  i-1  11 


a-1 


+  E  I(t.  exp{-  E  X?}  H  (1-hU  X?”1) 

j-1  J_J-  J  k-j+1  *  k-J-f-1  K  K 


,o-l 


(l-h»Jj(Xj+tj-1-T)“  *]  exp{-(X:j+tj-1-T)'*} 


0-1, 


+1{tn-l<TStn}  tl4aJn<Xn+tn-l-T)  1  «P{’(Xn+tn-rT)  } 

+I{t>t  }  , 

Q 

where  I(A}  is  the  indicator  function  of  the  set  A. 


(3-2) 


N 
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Generally,  let  £(t)  denote  the  prior  p.d.f.  of  r  .  Then,  the 
posterior  p.d.f.  of  x  ,  given  X^  ,  is 

S(T|X(n))  -  5(x)L(T;X(n))/D(5;X(n))  , 

where 

D(5;X(n))  -  /L(t;X(n))dP5(TSt)  . 

In  the  present  study  we  consider  a  prior  distribution  with  p.d.f. 


(3.3) 


(3.4) 


(l-p^e^-V 


,  if  TST. 


,  if  T>T. 


(3.5) 


where  0<p<l  and  0«j ><»  .  is  a  time  point  such  that,  with  high 

prior  probability,  1-p  ,  the  true  wear-out  point,  r  ,  exceeds  it.  We 
are  interested  therefore  in  time  points  greater  than  .  We  shall 
assume  accordingly,  without  loss  of  generality,  that  tq**0  .  For  the 
prior  p.d.f.  (3.5)  under  consideration,  the  function  (3.4)  assumes  the 
form  n 

D(p,*;X(n))  -  p  exp{-  I  X®}  n  [l+aJ.x"'1] 

1-1  1  i-1  1  1 


+  (1-p)  E  exp{-  E  XT)  n  [1-hxJ.XT  *] 
j-1  k-j+1  *  k-j+1  K  * 


[i-w,®,--)"-1)  «-<xrur^ 

0  3  J 

+  (l-p)*e  <l<tn-1  /V*U  fl+aJ  (X  -u)®'1]  e_(Vu)  du 
0  Q  n 

■*tn 

+  (l-p)e  n 

The  posterior  probability,  after  observing  ,  that  the  transition 

to  Phase  II  has  already  taken  place  (i.e.,  {tst^})  is  pQ  *  l-qQ  , 
where 


o-l,  -(X.-u)c 


(3.6) 


qQ  -  (l-p)e  n/D(p,4»;^n)) 


(3.7) 


If  pn  is  large  ve  have  high  evidence  that  the  change  has  already  taken 
place.  In  addition,  after  observing  the  Bayes  estimator  of  t  , 

for  squared  error  loss  function,  is  the  posterior  expectation  of  r  , 


given  jp  '  .  This  estimator  is  determined  by  the  formula 
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Tn(p,i|i;X(n*)  -  E{T|X(n\p,<J»}  - 

n-1  n  ,  n 

{ (l-p)i)»  E  n  [1-taJ.  Xf”1’]  exp{-  E  XT>  . 
j-1  k-j+1  K  k-j+1  K 

— il»t ,  ,  X.  .  ,  -(X.-u)a 

•e  j“1  /j  (t,  1+u)e"^,u  [1+aJ.  (X.-u)*'1]  e  j  du 

0  J”1  J  3 

+  (l-p)^/e  n_1  /n  (t  ,+u)e  [1+aJ  (X  -u)a  e  ^Xn  U^du 

q  n~i  n  n 

00  f  \ 

+  (l-p)i pe  Q  /  (tQ+u)e  ^udq}  *  D(p,i|i;X^n  )  .  (3.8) 

0  •  /_\ 

Notice  that  if  A*1  we  determine  PQ  and  Tn(p,t|);X'  ')  by  substi¬ 

tuting  in  the  above  formula  XX.  for  X.  and  i| >/\  for  ip  . 

A  A  X 

Thus,  we  can  compute  p  and  t  adaptively,  after  each  replace- 

HQ  A 

rnent,  and  decide  that  the  change  has  occurred  at  the  first  time  r  £t 
or  pQ  is  sufficiently  large.  Shiryaev  [  2 ]  studied  the  problem  of 
the  quickest  detection  of  the  occurrence  of  a  change  in  a  speclfided 
distribution.  He  showed  that,  when  the  distributions  before  and  after 
the  change  are  known,  the  optimal  Bayes  sequential  stopping  time  is  at 
the  smallest  nil  for  which  pQiII  .  H  is  a  threshold  that  can  be 
determined  to  minimize  the  total  expected  loss.  Our  problem  is  however 
different.  First,  we  do  not  have  just  two  known  distributions  with  a 
switch  from  one  to  another  at  an  unknown  change-point ,  t  ,  but  a 
sequence  of  wear-out  distributions  with  unknown  parameters  ,  and  we 
wish  to  detect  the  time  point  at  which,  for  the  first  time,  . 

Moreover,  we  do  not  have  a  sequential  stopping  problem,  but  an  adaptive 
decision  problem.  The  replacement  process  is  continuing  also  after  the 
decision  that  t^sO  has  been  reached.  The  implications  of  this  deci¬ 
sion  are  on  the  adaptation  of  the  scheduled  replacement  period,  A  , 
and  on  the  Inventory  of  spare  components.  In  the  following  section  we 
provide  some  recursive  formulae  for  the  determination  of  D(p,^;^^n^>  , 

A 

pn  *nd  Tn(p’*;V  * 
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4.  Recursive  Formulae 

From  (3.3)  and  (3.5)  one  immediately  obtains  that,  for  every 
t2tQ  (n-1,2,...) 


P{T2t|X^} 


(1-p)  e~  iflt 

D(p,^;X(n) 


(4.1) 


Confining  this  result  with  (3.7)  we  obtain  for  all  t2t 


P{T2t|Xw}  -  (l-pn)e 


(4.2) 


where  p  is  the  posterior  probability  of  {t£t  }  .  This  yields  that, 
d  /  .  n 

after  observing  X^  '  ,  the  posterior  distribution  of  t  for  tit^  is 
of  the  same  form  as  (3.1),  with  p  replaced  by  p  ,  i.e.,  the  poster¬ 
ior  p.d.f.  of  t  ,  given  satisfies 


(1-Pa)$e 


“Ip  (  T“t  ) 


,  if  xst 


,  if  T>t 


(4.3) 


According  to  (3.6),  the  function  =  D(p,tJi;X'‘  ')  ,  n2l  ,  can  be  com¬ 
puted  by  the  following  recursive  formula 


D  -  D  (1-KxJ  X®'1)  exp{-X®}  +  (l-p)e 
n  n-i  n  n  n 


O  e  u(l+aJ  (X  -u)a  ;  exp{-(X  -u)a}du 
0  n  n  n 


-  (l+aJX®'1)  exp{-X®}]  +  (l-p)e  n 


where  Dq  =  1 


Notice  that  for  all  X  ,  0<X<®  , 

D(p,*;X(n>)  -  D(p,*/X;XX(n)) 


(4.4) 


(4.5) 


Employing  (3.7)  and  (4.4)  we  obtain 


n  'k  n-1  D  .  -*X 

U~P)e  .  .  n-i  en 

D  -  D 

Qrl  n 


^1-Pn-l^  D  Id  T 

n  n-1 


(4.6) 
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*2<p,*,x)  -  pe  *  (l+2Jx)  + 

(l-p)*«“*x*l|,2/4  ft  [*(^(x  -  4-  )>  +  *(  )  -1]  (1-M*) 

2  ft 

+j(i-P)^e-^2/4  /nr  [♦(  -t- )  -  ♦<  -JL  -  /2x>] 

ft.  ft 


+(l-p)e 

where  $(z)  and  4(z)  are,  respectively,  the  p.d.f.  and  the  c.d.f.  of 
the  standard  normal  distribution. 

We  discuss  now  the  recursive  determination  of  the  Bayes  estimator 
of  t  .  According  to  (3.8) 


(4.12) 


(4.13) 


where 


r  n-l  n  .  n 

E  -  1  (1-p)*  E  n  [14«J  x“_i]  exp{-  Z  X?) 

1  j-1  k-j+1  k  k  k-j+1  * 

-*t  _  Xj 

e  ^  1  /  (tj_j+u)e  ^u[l-hxJj  (Xj-u)a_1]  exp{-(Xj-u)a}du 


u  X 
.  n 


+  (l-p)^e  n  ^  /  (t  -+u)e  *^U  [1+aJ  (X  -u)°  *]  exp{-(X  -u)a}dul  *  D 
0  n-l  n  n  n  J  u 


Define  the  function 


n  (p,«|»,t,x)  -  (l-p)e”’*%  /  (t-y)  (l+aJya~1)e^y~y  dy 

0 


Then, 


Ei  “  na(p,'l>,t1,x1)/D1 


(4.14) 


(4.15) 


(4.16) 


and,  for  every  n2l 


s«u ' 


♦  *  D»+i 


(4.17) 
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The  function  ^(Pf <)'» t,x)  can  also  be  written  in  terms  of  the 
.d.f.  and  c.d.f.  of  the  standard  normal  distribution  as 

n2(P»t>t>x)  *  (l-p)<|>e  *t+*  ^ 

2 

{(4(/f(x  -  4-  ))+*<-£->-  u  (*: - 2~  ^  J(t+ - 2 - 1)) 

L  /l 

+  4>(  — )  (tJ/I - — - J  "•(r 1  )  —  4>  C  ~ 

/2  /2  /2  /2  /2 

-  J(  — +  /2x))} 

/2 
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5.  Some  Numerical  Simulations 

We  present  here  results  of  some  simulations,  which  Illustrate  numer¬ 
ically  the  characteristics  of  the  procedure  developed  In  the  previous 
sections.  In  order  to  simulate  a  random  variable  Y  ,  having  a  wear-out 
distribution  with  parameters  A,  a  and  t  ,  proceed  as  follows.  First, 
simulate  a  random  variable  U  ,  having  a  uniform  distribution  on  (0,1)  . 
Then,  solve  for  Y  the  equation 

exp  {-AY  -  Aa(Y-r)“}  -  U  .  (5.1) 

Finally,  determine  X  -  min  (Y,A) 

In  the  present  study  we  restrict  attention  to  the  case  of  a*2  . 

Accordingly , 

-  j  In  0  ,  if  U*e~XT 

Y  -  (5.2) 

T‘ir+x(i“lnU'  » if  D<e”XT 

In  the  following  examples  we  consider  a  system  with  -  200  [hr]  , 

A  «  225  [hr]  ,  t  ■  750  [hr]  ,  a  ■  2  and  the  prior  probability  that 

{tsO}  is  p  ■  .2  .  For  the  simulation  of  the  random  variable 

X  ,  (n*l,2,...)  we  apply  a  value  of  t  -  (t-t  .) ,  .  The  values  of 
n  n  n-i.  + 

the  posterior  probabilities,  pQ  ,  and  those  of  the  Bayes  estimators, 

,  are  computed  adaptively  after  each  stage.  In  Table  5.1  we  present 
the  results  of  such  a  simulation  run. 


Table  5.1.  A  Simulation  Run  with  Parameters 
A  -  1/200,  ift  -  1/1000,  p  -  .2,  t  -  750 


n 

K 

n 

J 

n 

t 

n 

pn 

a 

Tn 

1 

41.49 

1 

41. 

0.288 

742.71 

2 

156.56 

1 

198. 

0.473 

646.85 

3 

108.39 

1 

306. 

0.631 

514.63 

4 

225.00 

0 

531. 

0.437 

921.39 

5 

225.00 

0 

756. 

0.309 

1317.66 

6 

31.19 

1 

788. 

0.385 

1237.91 

7 

23.03 

1 

811. 

0.446 

1170.46 

8 

115.07 

1 

926. 

0.610 

1011.46 

9 

114.04 

1 

1040. 

0.742 

866.27 

10 

20.53 

1 

1060. 

0.779 

815.35 

11 

27.09 

1 

1087. 

0.820 

758.42 
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We  see  In  the  above  example  that  the  first  value  of  t  larger  than 

x  is  t  ■  1040  [hr]  .  This  is  greater  than  t  by  290  [hr]  ,  which 
n  n 

is  a  little  over  than  1  MTBF.  The  corresponding  value  of  p  is 
*  »  ® 

p  *  .742  .  The  estimate  of  t  is  t  ■  866.27  .  We  see  also  that  if 
n  n  * 

we  defer  the  decision  until  pQ£.8  then  t  -  t  ■  337  [hr]  and 

t  ■  758.42  .  Thus,  the  second  decision  (stopping)  time  provided  a 
n 

more  accurate  estimate  of  r  .  In  Table  5.2  we  present  frequency  dis- 

*  * 

tributions  and  the  exact  means  and  standard  deviations  of  W  ■  t  -  x 

n  n 

as  obtained  by  M"100  independent  simulation  runs,  for  each  one  of  the 

prior  parameters  <j»  -  1/500,  1/750  and  1/1000  .  The  decision  time  in 

these  simulations  is  the  first  value  of  t  greater  than  t  ,  i.e., 

*  n  n 

t  »  least  n^l,  such  that  t  2x  (5.3) 

n  n  n 


Table  5.2.  Frequence  Distributions,  Means  and  Standard 

* 

Deviations  of  W  in  M*100  Simulation  Runs 
-  n  —  - ■ - - - 


Mid  Point 
Interval 

Frequencies 

l/ti(-50  0 

l/it-750 

1U-1000 

-500 

10 

0 

0 

-400 

35 

28 

15 

-300 

17 

37 

22 

-200 

21 

7 

19 

-100 

8 

11 

6 

0 

0 

0 

10 

100 

6 

7 

10 

200 

3 

9 

5 

300 

0 

1 

7 

400 

0 

0 

6 

Mean 

St.  Dev. 

-278.15 

167.44 

-218.44 

201.88 

-11 '■'.72 

239.82 

* 

We  see  in  Table  5.2  that  the  decision  time  tQ  tends  to  yield  too  many  early 
decision  points.  Zf  l/ip»l,000  [hr]  the  results  are  significantly 
better  than  in  the  case  of  l/ij>“500  or  750  [hr]  . 

The  situation  seems  to  be  better  in  the  case  of  the  decision  times 
defined  by 


t  ”  least  n  s.t.  p  i.9 
n  n 

**  - 

In  Table  5.3  we  present  the  frequency  distributions  of  Wq 
obtained  In  independent  simulation  runs. 

Table  5.3.  Frequency  Distributions.  Means  and  Standard 

M 

Deviations  of  V  in  M*1Q0  Simulation  Runs 


Mid  Point 
of 

Intervals 


Frequencies 

1/^500 

1/^750 

1/^-1000 

Mean 
St.  Dev. 


65.60 

311.78 


152.85 

325.40 


202.13 

322.37 


It  seems  that  the  Bayes  decision  times  tQ  based  on  the  prior 
distribution  (3.5),  with  a  proper  choice  of  the  <|*  value  (not 
too  small) ,  provide  good  results  in  the  problem  of  detecting  a 
shift  to  the  wear-out  phase. 
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